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Abstract 

We study the synchronization of nearest neighbors coupled oscillators in a ring. We derive an 
analytic form for the phase difference among neighboring oscillators which shows the dependency 
on the periodic boundary conditions. At synchronization, we find two distinct quantities which 
characterize four of the oscillators, two pairs of nearest neighbors, which are at the border of the 
clusters before total synchronization occurs. These oscillators are responsible for the saddle node 
bifurcation, of which only two of them have a phase-lock of phase difference equals ib7r/2. Using 
these properties we build a technique based on geometric properties and numerical observations 
to arrive to an exact analytic expression for the coupling strength at full synchronization and 
determine the two oscillators that have a phase- lock condition of ±7r/2. 
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I. INTRODUCTION 



Coupled oscillators have been used to understand the behavior of systems in physics, 
chemistry, biology, neurology as well as other disciplines. In particular, they are used to 



timode lasers, vortex dynamics 
3|]. These systems have been 



model phenomena such as: Josephson junction arrays, mu 
in fluids, biological information processes, neurodynamics 
observed to synchronize themselves to a common frequency, when the coupling strength 
between the oscillators is increased Is-Sj. Although all these phenomena have different 
dynamics, their synchronization features might be described using a simple model of weakly 
coupled phase oscillators such as the Kuramoto model and variations to adapt it for finite 
range interactions which are more realistic to represent physical systems sj-Gj. But, finite 
range interactions difficult the analysis and search for analytical solutions. In spite of that, 
in order to figure out the collective phenomena when finite range interactions are considered, 
it is of fundamental importance to study and to understand the case of nearest neighbor 
interactions, which is the simplest form of the local interactions. In this context, a simplified 
version of the Kuramoto model with nearest neighbor coupling in a ring topology, which we 
shall refer to as locally coupled Kuramoto model (LCKM), is a good candidate to describe the 
behavior of coupled systems with local interactions. The LCKM has been used to represent 
the dynamics of a variety of systems. Specifically, it has been shown that a ladder array 
of Josephson junctions can be expressed by a LCKM [7]. Phase synchronization, in nearest 
neighbors coupled Rossler oscillators and locally coupled lasers where local interactions are 
dominant, can also be described by the LCKM jsHl. Other examples are: the occurrence 
of travelling waves in neurons, chains in disorders, multi cellular systems in biology, the 
dynamics of an edge dislocation in a 2D lattice and an antenna array in communication 
systems 

While in the Kuramoto model of long range interactions one has to get a solution in 
a mean field approximation, in the local model it is necessary to study the behavior of 
individual oscillators in order to understand the collective dynamics. Due to the difficulty 
in applying standard techniques of statistical mechanics in order to obtain a close picture 
of the effect of the local interactions on synchronization, we rely on a simple approach 
to understand the coupled system with local interactions, by means of numerical study of 
the temporal behavior of the individual oscillators. In this case, numerical investigations 
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provide a good tool to understand the mechanism of interactions at the stage of complete 
synchronization which in turn helps to get an analytic solution. Earlier studies on the LCKM 
show several interesting features including tree structures with synchronized clusters, phase 



slips, bursting behavior, saddle node bifurcation and so on |15l. |16|. It has also been shown 
that neighboring elements share dominating frequencies in their time spectra, and that this 
'eature plays an important role in the dynamics of formation of clusters in the local model 



13, 



18| : that the order parameter, which measures the evolution of the phases of the nearest 



neighbor oscillators, becomes maximum at the partial synchronization points inside the tree 



of synchronization 



19| and a scheme has been developed based on the method of Lagrange 



multipliers to estimate the critical coupling strength for complete synchronization in the local 
Kuramoto model with different boundary conditions [20]. In addition, based on numerical 
investigations, we identified two oscillators which are responsible for dragging the system 
into full synchronization 2l|, and the difference in phase for this pair is ±7r/2. These two 
oscillators are among two pairs of oscillators which are formed by the four oscillators at the 
borders between major clusters in the vicinity of the critical coupling. Using these findings 
we developed a method to obtain a mathematical expression for the estimated value of the 
critical coupling at which full synchronization occurs, once a set of initial conditions for the 
frequencies of the N oscillators is assigned [22 1. 

In this work, we use the fact that at the stage of full synchronization, all oscillators have 
a common frequency and the time dependence of the phase difference among neighboring 
oscillators will vanish. Using this we are able to derive an analytic expression for the phases 
of oscillators, which shows the actual dependence on the periodic boundary condition. It is 
clear that this effect will decrease when the number of oscillators N increases, and the result 
will converge to that of a free chain as tends to infinity {22!. Even so, both problems 
are different and have interesting peculiarities when the number of oscillators is finite. This 
problem is not a mere detail in a theoretical problem since a ring with a finite number 



of osci 
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23 



lators has many applications in electronics, coupled lasers and in communications 



24| . In the process of finding the solution, using the saddle node bifurcation which 



dominates the dynamics at synchronization, we come across two quantities which will permit 
us to identify the four oscillators at the borders between major clusters, where only two of 
them will have a phase-lock condition of |vr/2|. We use the properties of these oscillators 
and the ones at the boundaries to define a triangle and use its trigonometric properties to 
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derive an analytic formula for the critical coupling. Finally, we can identify directly the pair 
of oscillators which has the phase-lock condition, which depends only on the set of initial 
frequencies. 

This paper is organized as follows: In Section II, we introduce the LCKM with boundary 
conditions. We determine analytically the critical coupling at the stage of complete syn- 
chronization. Finally, in Section III we give a conclusion which is based on a summary of 
the results. 

II. THE MODEL 

The LCKM can be considered as a diffusive version of the Kuramoto model, and it is 

expressed as 

9i^0Ji + K (sin (pi - sin (1) 

with periodic boundary conditions 9i+iy = 6i and phase difference (pi = ^^j+i — 9i for i — 
1,2, ...,N. The set of the initial values of frequencies Ui arc the natural frequencies which 
are taken from a Gaussian distribution and K is the coupling strength. These nonidentical 
oscillators of system 1 cluster in time averaged frequency, under the influence of the coupling, 
until they completely synchronize to a common value given by the average frequency ouo — 
Jf X^ili ^ critical coupling as shown in Figure 1. At the vicinity of Kc, there are only 
two clusters of successive oscillators whose borders are nearest neighbors. The oscillators 
remain synchronized for K > Kc when all phase differences and frequencies are constants. In 
Figure 1, we show the synchronization tree for a system with iV = 40 oscillators, where the 
elements which compose each one of the major clusters, just before complete synchronization 
occurs, are indicated in each branch. These clusters merge into one at Kc where all oscillators 
have the same frequency. The major clusters, at the onset of synchronization, just before 
Kc contain Ni and N2 oscillators, where N — Ni + N2. It is not necessary for these clusters 
to have the same numbers of oscillators. We point to the oscillators at the border of each 
cluster, which are formed with successive elements, thus their bordering elements are nearest 
neighbors, as mentioned above; i.e, + 1,171 and m + l. An interesting fact emerges here: 
there is a phase- locked solution, where the phase difference between two oscillators is ±7r/2, 
and it is always valid for one and only one phase difference. This can be the difference 



4 



between phases of any two nearest neighbors of the four oscillators at the border of the 



clusters at the onset of synchronization 



2l|. However, we can not directly allocate these 



two oscillators unless we do it numerically. At this point, we remind the reader that the 
location of each oscillator corresponds to a well defined entity, characterized by an initial 
frequency, therefore wondering about the meaning of the location of these borders as well 
as the position of the boundaries, corresponds to knowing which of the many oscillators in 
our systems will dominate the dynamics. 

Thus, the time evolution of the phase differences among neighboring oscillators will be 
written as 

^i = Ai- 2Ksm{(f)i) + irsin(0,_i) + sin(0i+i), (2) 

where Aj = Wj+i — coi. We use the fact that at Kc, the quantities (pi = 0, to derive the 
following expression 

sin(0i) = -rr + sin(07v), (3) 



where, if, 



(N-i+l) 
N 



According to equation 3, we need to 



identify not only the two oscillators, say j and j + 1, which will phase-lock with (pj = |vr/2|, 
but also the value of the phase difference between the oscillators identified as boundaries to 
find the critical coupling. In order to tackle this difficulty we depend on numerical simula- 
tions of system 1 for different number of oscillators and sets of initial frequencies Ui, as well 
as on studying the quantity Hi. It should be noted that the quantity max {|-?ij|} determines 
the value of the critical coupling and the oscillators locked on |7r/2| at synchronization, for 
the case of a chain of free ends (open boundaries). This can be clearly understood from 
equation 3, when sin((/)Ar) = 0, and | sm{(f)j)\ = 1, which coincides with the same index of 
max{|ifj|}, and henceforward = max{|ifj|}. However, when periodic boundary condi- 
tions are turned on, this fact does not hold true, since sin(0Ar) can be either positive or 
negative. Therefore, its presence in equation 3 modifies the quantity Hj which matches 
the same index of the phase difference which locks to |7r/2|, and the oscillators no longer 
synchronize to a single frequency Uo at a critical coupling defined as max{|/ij|}. Helped by 
numerical simulations, we can find the four oscillators at the edges of the two clusters at 
the onset of synchronization, in the vicinity of Kc, and see that only one phase difference is 
locked to |±7r/2| at as well as to determine which two oscillators among these four are 
responsible for the saddle node bifurcation. 
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When we study numerically the phase differences among neighboring oscillators we find 
two distinct values among all phase differences, which are characterized by being the absolute 
maximum and minimum among all phase differences and we shall see that they correspond 
to the two phase differences of the four oscillators at the borders [22]. But, only one of them 
will have the phase- lock condition of | ±7r/2| at Kf.. If we now proceed to study the behavior 
of the quantity Hi, we see that the maximum and minimum values of Hi always coincide 
with the same indexes of the maximum and minimum of sin({/)j), respectively. Therefore, 
based on the numerical simulations, we can distinguish two quantities Hi and Hm and their 
corresponding phase differences 0^ and (f)m- Any one of them, H^ or H^ can be the absolute 
maximum or the absolute minimum. This holds true for the two values sin(0£) and sin(0m) 
also. Comparing between the absolute values of either the two quantities H^ and H^ or the 
two quantities sin(0£) and sin(0m), we can not say which oscillators will have the phase-lock 
condition. However, we find always, and for any with any sets of Wj, that the signs of Hg^ 
and sin(0^), as well as the corresponding Hm and sin(0m) are always the same, independent 
of being positive or negative. Figure 2a shows the plot of selected values of sin(0j) versus 
time at for the same realization of Figure 1. We can see that sin(0£) = 1 for £ = 7 and 
sin(0m) < 0, for m = 23, which is the absolute minimum among all corresponding values 
of all oscillators. Figure 2b shows the values of Hi versus the positions of oscillators in the 
chain for the same case of Figure 1, where we clearly see that H^ > for £ = 7 is the absolute 
maximum, while Hm < for m = 23 is the absolute minimum. Since these oscillators are 
obviously special, we shall focus on them in search for the exact solution of the problem. 
We see that the maximum and minimum of Hi determine two major cases: H^ > and 
Hm < 0. Once we select the first case, we notice that, again, we encounter two cases. We 
shall see that only one of them decides the exact value of the critical coupling. These two 
sub-cases are a) sin(0^) = 1 with sin(0m) > —1, and b) sin(0^) < 1 with sin(0m) = — 1- 
For case (a) we use three equations from system ([3]) for i,i — 1 and m, from which we can 
write '^'"^"^^•^ = = This relation correspond to the properties between the 

angles and sides of a triangle, where the angles and sides are defined as: a = —(pN + 7r/2, 
/3 = -0^ + 7r/2 and 7 = -0£_i + 7r/2, a = {HiY'\ b = {Hi-Hmf^ and c = {Hi-H,_,f'\ 
as shown in Figure 4. After some manipulation, we get an expression for (pm as 
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TABLE I: Values of K,, from simulation of system 1 and from equation 5. 



N 


Kc'. simulation 


Kc- equation ([5]) 


30 


3.73094125 


3.72862539 


40 


4.29473534 


4.28553971 


50 


4.48415639 


4.47935214 


100 


5.86827841 


5.86639415 


200 


7.96802973 


7.96457428 



- 2 cos^^ <^ \ . (4) 



"2 t 2ac 

Thus, the value of the critical coupling becomes: same method 

to case (b) (sin(0/) < 1 with sin(0m) = —I), we get i^"™ = -if^(^y- Therefore, the value 
of the critical coupling when Hi > and Hm < is 

ir, = max{irf,ir™}. (5) 

Following the same method for the case Hf < and Hm > 0, we find: K^. = and 
-^c" ~ • "^^^ value of the critical coupling is again given by equation 5. Table 1 

shows the results from numerical simulations of system 1 which are in good agreement with 
the values obtained from equation 5 for the same sets of initial frequencies. 

Equation 5 allows us to determine whether (pi or (pm has the phase- lock condition \ti/2\ 
at depending whether the selected value for Kc is Kc or K"^, respectively. After this is 
done we can determine the value of sin(0Ar) ^ 0, thus showing how the phase difference (pM, 
consequence of the periodic boundary conditions, influences the dynamics of the system. 
We notice from Figure 2a that the value of sin(04o) 7^ 0, and it modifies the result of the 
critical coupling, as expected from equation 3. The numerical simulations, for different N 
and for different sets of Wj, show that the system with periodic boundary conditions has a 
critical value which depends on A^, and that the influence of sin(07v) decreases as increases 
which may give indications that it becomes ineffective, but in this case the system never 
synchronizes at the critical strength given only by the absolute maximum of Hi, as it should 
be for a free chain. We find also, that sin(0Ar) remains different from zero up to A^ ~ 2000. 
Even though, sin(0iv) has a small value, but since it exists, it enforces the oscillators to 
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synchronize approximately as an average of Hi and Hm, which are not in generally equal to 
each other. For larger values of A^, we expect that \H£\ ^ \Hm\ and sm{(f)N) = 0, and hence 
the critical coupling is determined by the absolute value of Hi. As N ^ oo, the system of 
oscillators in a ring synchronizes at the same value of the coupling strength as a chain of 
free ends, as expected |15 |. 

Qualitatively, we can compare the cases of synchronization for a ring and for a free chain 
for finite number of oscillators in order to understand why the critical coupling depends 
on both values max {Hi} and min {Hi} for the first case while the free chain only depends 
on max{|i7j|}. If we consider both systems above full synchronization and decrease the 
coupling strength we observe that a saddle node bifurcation occurs at the critical value and 
they split into two clusters of unequally number of oscillators in general. The chain splits 
into two clusters as shown in figure 5a, where only two oscillators seem to have a new role 
in the dynamics, j and j + 1 are located in different clusters. This splitting occurs when the 
phase difference (pj = | vr/2 1 , and it is the index that maximizes max { | -ffi | } which determines 
the value of Kc for the chain. The ring at the critical coupling splits into two clusters of 
unequal numbers of oscillators, as shown in figure 5b. The difference is that for the ring two 
sets of oscillators acquire a special role in the dynamics: i, £+1, m and m + 1 but only one 
phase difference will be 0^ = |vr/2| and it can be either 0£ or 0^, where the two indexes i and 
m match max {H^} and min {Hi}, respectively. A priori we do not know which will be. From 
the previous discussion we can write Kc = \Hj ± where j = i or m and 6 depends 

on the periodic boundary conditions, in order to get |sin(0j)| = |Xj • + sin(0Ar)| = 1 and 
Xj = Hj/Kc 7^ ±1 and hence sin(0Ar) "decides" which will be the pair of oscillators which 
phases-lock at |'/r/2|. Therefore, for finite A^, the oscillators synchronize at a critical coupling 
that depends on both quantities max {Hi} and min {Hi}, and this is due to the presence of 
the periodic boundary conditions. From equations 3 and 5 we can write: Kc = j^y^^^j^^; 
with I sin(0j)| < 1, where j = m and | sin(0£)| = 1, and j = i when | sin(0m)| = 1- 



III. CONCLUSIONS 



We have studied the synchronization of coupled oscillators in a locally coupled Kuramoto 
model with nearest neighbors coupling with periodic boundary conditions. Particularly, we 
analyzed the system of oscillators at the stage of complete synchronization. We see that the 
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dynamics is determined by a particular set of oscillators which are located at the borders 
of the major clusters which will meet at the critical coupling to form one cluster of all 
synchronized oscillators. Using the boundary conditions we can make a correspondence of 
the values of some definite quantities Hi with the sides of a triangle with the aid of some 
trigonometric properties. These quantities correspond to those of the border oscillators and 
their nearest neighbors. Using these two quantities, we derive a mathematical expression 
for the critical coupling at synchronization. In addition, we are also able to determine the 
two oscillators which will have a phase- lock condition of |7r/2|. All these properties can be 
calculated a priori since they depend only on the set of initial frequencies Wj. 
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FIG. 1: Synchronization tree for a system of 40 oscillators with detailed composition of each cluster 
before full synchronization. Here we point the oscillators at the border 
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FIG. 2: Selected values of sin(^i) at Kc for the oscillators of Figure 1. 
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FIG. 3: The values of Hi versus oscillator position, for the oscillators of Figure 1. 
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FIG. 4: Triangle of sides: a = (He)^/^, b = {He - Hmf'^ and c={Hi- Ht^if^. 
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Fig. 5(b) 



FIG. 5: The two major clusters at the critical coupling in (a) the chain of free ends and (b) the 
ring. 



16 



